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ABSTRACT 


In  order  to  obtain  the  cross-scan  position  of  an  optical 
targe!;,  more  than  one  scanning  detectors  are  used.  As  expected, 
the  cro&s-scan  position  estimation  performance  degrades  when  two 
nearby  optical  targets  interfere  with  each  other.  Theoretical 
bounds  on  the  two-dimensional  parameter  estimation  performance 
for  two  closely  spaced  optical  targets  are  found.  Two  particular 
classes  of  scanning  detector  arrays,  namely,  the  crew's  foot  and 
the  brickwall  (or  mosaic)  patterns,  are  considered. 
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INTRODUCTION 


Recently,  considerable  attention  has  been  focused  upon  the 
accuracy  of  parameter  estimates  of  an  optical  system  for  closely 
spaced  objects  (CSOs)  [1]  -  [7].  The  intensity  and  location  in¬ 
formation  of  an  optical  image  on  the  sensor  focal  plane  is  collect¬ 
ed  by  a  scanning  detector  as  shown  in  Fig.  1.  The  detector  output, 
a  time  function,  is  a  convolution  of  the  optical  image  (which 
contains  intensity  and  location  information)  and  the  detector 
response.  There  are  various  types  of  noise  sources  at  the  output 
of  the  detector,  therefore  estimating  the  image  position  and  in¬ 
tensity  from  the  noisy  output  can  not  be  errorless.  Theoretical 
lower  bounds  on  the  accuracy  of  parameter  estimates  have  been 
derived  in  [1]  -  [4], 

The  Cramer- Rao  lower  bound  technique  is  used  in  [2]  -  [4]  to 
calculate  the  best  achievable  performance  of  any  unbaised  estimator 
of  the  unknown  parameters.  The  results  presented  in  [1]  are 
obtained  using  a  different  error  analysis  technique.  Other  than 
some  numerical  problems  pointed  out  in  [4],  which  occurred  at 
small  target  separations,  the  approach  presented  in  [1]  should 
have  results  identical  to  those  obtained  by  the  Cramer-Rao  lower 
bound  approach.  There  are  algorithms  developed  to  perform  the 
parameter  estimation  for  closely  spaced  objects  from  noisy 
measurements,  for  example  [5]  -  [7].  Monte  Carlo  simulation 
results  are  presented  in  (5]  and  [7].  In  [5],  the  well  known 
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Fig.  1.  Closely  spaced  object  (CSO)  resolution  problem  for 
optical  system  with  scanning  detector. 


maximum-likelihood  estimator  is  implemented  for  parameter  estimation 
and  the  Akaike  information  criterion  [8]  is  employed  to  estimate 
the  number  of  targets  in  the  CSO  cluster. 

All  analytical  and  simulation  results  presented  in  the  past 
pertain  only  to  the  in-scan  direction  and  assume  a  vertical 
rectangular  detector  as  shown  in  Fig.  1.  Many  detector  array 
arrangements  have  been  considered  by  sensor  designers  to  obtain 
position  information  of  a  single  optical  target  in  the  cross-scan 
direction.  Some  typical  detector  patterns  considered  are  shown  in 
Fig.  2.  The  time  between  detector  crossings  for  the  chevron  and 
crows 's  foot  patterns  can  be  used  to  determine  the  cross-scan 
position  of  the  target.  With  two  targets,  this  problem  becomes 
more  difficult.  Figure  3  shows  three  individual  detector  outputs 
of  a  crow's  foot  pattern  for  different  orientations  of  target 
locations.  It  is  obvious  that  with  one  linear  detector  (detector 
M)  output  alone,  one  can  not  obtain  cross-scan  position  information. 
From  the  results  we  have  presented  in  [2]  -  [5],  one  can  obtain  in¬ 
scan  position  information  of  two  CSOs  no  matter  how  close  they  are 
if  the  detector  output  is  noise  free.  It  is,  however,  not  true  for 
a  noisy  output;  the  estimation  errors  increase  drastically  as  the 
CSO  separation  decreases.  For  the  chevron  pattern  (detectors  L 
and  R) ,  the  cross-scan  position  information  can  be  obtained  with¬ 
out  error  from  noise  free  outputs  for  all  target  locations.  For 
noisy  outputs,  there  are  target  locations,  for  example,  case  II  of 


Typical  detector  patterns  have  been  considered 
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Fig.  3.  Individual  detector  outputs  of  a  crow's  foot  pattern 
for  different  orientations  of  target  locations. 
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Fig.  3,  where  the  target  position  estimation  errors  may  be  very 
large  no  matter  how  large  the  actual  target  separation  is.  This 
is  because  there  is  a  CSO  situation  for  the  detector  R.  A  more 
complicated  detector  array  arrangement  may  solve  this  problem, 
for  example  the  crow's  foot  pattern  shown  in  Fig.  3.  The 
achievable  estimation  performance  of  the  cross-scan  position 
of  two  CSOs  for  these  particular  detector  patterns  has  not  been 
studied  analytically  previously.  It  is  the  purpose  of  this  report 
to  derive  the  Cramer-Rao  lower  bounds  on  the  variance  of  parameter 
estimates  as  a  function  of  both  scan  and  cross-scan  separations 
for  CSOs  measured  by  detector  patterns  shown  in  Fig.  2.  The 
performance  bounds  presented  in  this  report  are  derived  with  the 
assumption  that  all  detector  outputs  are  utilized  optimally  for 
unbiased  parameter  estimation;  no  particular  signal  processing 
procedure  has  been  adopted  in  the  analysis. 

The  general  formulation  of  this  problem  and  its  associated 
Cramer-Rao  bounds  are  presented  in  the  next  section.  In  Section 
III,  results  for  linear,  chevron,  and  crow's  foot  arrays  using  a 
particular  optical  point  spread  function  are  presented.  For  other 
point  spread  functions  and  detector  widths,  one  can  utilize 
the  formulas  provided  in  Appendix  A  to  calculate  the  associated 
estimation  performance  bounds.  In  Section  IV,  the  brickwall  (or 
mosaic)  detector  patterns  are  considered.  Some  typical  results  foi 
the  same  point  spread  function  considered  in  Section  III  are  pre- 
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sented.  Formulas  that  can  be  used  to  calculate  the  estimation 
performance  bounds  are  provided  in  Appendix  B.  A  summary  and 
conclusions  are  given  in  Section  V. 

II.  ANALYSIS 

Let  us  assume  there  is  an  array  of  M  detectors  in  the  focal 
plane  of  a  scanning  optical  sensor;  some  typical  patterns  of  these 
detectors  are  shown  in  Fig.  2.  Let  (x^y^)  he  the  center  and  0^ 
be  the  orientation  of  the  kth  detector  with  the  center  of  the 
focal  plane  scanning  along  the  x-axis  as  shown  in  Fig.  4.  Let 
Pk(t,x,y)  be  the  output  of  the  kth  detector  to  a  unit  strength 
point  source  located  at  (x,y)  at  t=ot  A  point  source  is  called  a 
unit  strength  point  source,  if  for  (x,y)=(xk,yk)  and  6k=90°,  the 
output  signal  Pk(t,x,y)  satisfies: 

/QO 

Pk2 tt,x,y)dt  -  1.  (1) 

CO 

it  is  important  to  note  that,  Pk(t,x,y)  is  a  convolution  of  the 

point  source  blurred  image  (e.g.,  point  spread  function,  PSF,  ror 

a  monochromatic  source)  on  the  focal  plane  and  the  response  of  the 

kth  detector.  In  general,  it  is  a  complicated  function  of  the 

detector  location  (xk,  yk)  shape,  orientation  ( ©k ) ,  and  the 

instantaneous  target  location  (x-t,y)  on  the  focal  plane  at  time  t. 

A  single  point  source  located  at  (x,y)  with  amplitude  a  will 

tin  the  remainder  of  this  note,  a  point  source  (or,  target)  location 
is  always  referred  to  its  location  on  the  focal  plane  at  t=0. 


Fig.  4.  Individual  detector  location  with  respect  to  the  center 
of  focal  plane  and  its  orientation  with  respect  to  the  scanning 
direction. 


have  a  noiseless  output 


sk(t)  =  apk(t,x,y)  (2) 

at  the  kth  detector  at  time  t. 

The  problem  treated  here  concerns  the  measurements  of  a  pair 
of  point  sources  from  the  output  of  an  M  detector  array.  Let 
( a (xlf  y j ) )  and  (a2,  (x2,  y2))  be  and  location  of 

the  point  sources.  The  noiseless  output  of  the  kth  detector  can 
be  written  as  the  following: 

sk(t)  =  a1pk(t,x1,y1)  +  a2pk (t,x2,y2) .  (3) 

We  wish  to  determine  (a^,  (x^/y^))/  i=l,2  from  noisy  measurements, 
taken  at  the  output  of  each  detector,  i.e., 

yk(t)  =  sk(t)  +  k=l,2,...,M,  (4) 

where  nk(t)  is  a  white  Gaussian  noise  with  two-sided  power  spectrum 
density  Nq/2.  Furthermore,  these  noises  are  mutually  independent, 
that  is 

E  [n^ (t ) n j (t ) 3  =  0,  for  all  t  and  t,  when  i^j.  (5) 

The  Cramer-Rao  lower  bound  on  the  variance  of  any  unbaised 
estimator  of  the  unknown  parameters:  amplitude  of  object  1,  a^ , 
amplitude  ratio,  R  =  a2/a^,  location  of  object  1,  (x1,y^) ,  separa¬ 
tion,  r,  and  orientation  of  objects,  0  (the  angle  between  the 
detector  scan  direction  and  the  vector  (x2,y2)  -  (5^,^),  as  shown 
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in  Fig.  5,  can  be  obtained  by  inverting  the  Fisher  information 
matrix,  F,  whose  (i,j)th  element  is  [9] 


=  E 


3  In  A 


9  In  A| 

3  a.  } 


(6) 


where  ou  denotes  the  ith  unknown  parameter,  namely,  oc^a^,  ct2=:R' 
a3=*lr  a4=^l'  a5=r  anc*  a6~0/  an<^  w^ere  ^  is  t^ie  likelihood 
ratio  [9,  p.  274] 


and 


In  A 


s(t)  = 


/•  00 

f  y1 (t) s (t) dt  ■ 

_oo 

-/=T< 

°( * 

•  00 

•'  —  oo 

"s^t  r 

s2(t) 

• 

,  y(t)= 

y2(t) 

• 

• 

-3M(tL 

• 

JM(tL 

(7) 


(8) 


Substituting  (7)  and  (8)  into  (6),  and  using  the  statistical  model 
of  the  problem  described  in  (4)  and  (5) ,  we  have 


Having  an  explicit  expression  for  the  output  function  pk (t,xi,yi) , 
i=l', 2  and  k=l,2,...,M,  one  can  compute  F^  according  to  (9).  In 
the  following  sections  we  will  consider  two  simple  cases  where  the 
partial  derivatives,  3  sk(t)/3  ou,  i=l,...,6,  can  be  obtained  explicitly. 
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III.  LINEAR,  CHEVRON,  AND  CROW'S  FOOT  DETECTOR  ARRAYS 


In  this  section#  a  special  class  of  detector  patterns  is 
considered.  The  detector  array  consists  of  M  parallogram  detectors 
with  the  shorter  width  along  the  scan  direction  as  shown  in  Fig.  2. 
The  detector  pattern  is  called  a  linear  array  if  all  0^*90° , 
a  chevron  array  if  6^  takes  only  two  values,  one  greater  than  and 
one  less  than  90°,  and  crow's  foot  if  @k  takes  three  values  (i.e., 
a  combination  of  linear  and  chevron  arrays) .  For  simplicity  in 
the  analysis,  let  us  make  the  following  assumptions s 

1)  the  center  of  each  detector  is  located  along  the  scan 
axis,  i.e.,  yk=0  for  all 

2)  the  output  of  each  detector  to  a  unit  strength  point 

4* 

source  located  at  (x,y)  has  the  same  functional  forms 

Pk(t,x,y)  »  p(t-(x-xk)  +  y  cot0k) ,  (10) 

where  p(*)  is  a  symmetric,  non-negative-valued  function.  Further¬ 
more,  let  p(t)  be  the  autocorrelation  function  of  p(t),i.e.. 


P  (t) 


p(T)p(x-t)dx. 


(11) 


It  is  easy  to  see  that  P ( 0 )  =1  from  the  definition  of  Pk(t,x,y)  in 

(1). 

Rewriting  (3)  in  terms  of  the  parameters  considered  in  Section 


tThis  implicitly  assumes  that  the  detector  has  infinite  length  or 
no  edge  effects  on  the  detector  output. 
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II  with  the  above  assumptions ,  we  obtain 

sk(t)  =  a1p|t-(x1-xJc)  +  yjcotej  + 

a^^Rp  |t- (Xj+rcosG-x^)  +  (y-^+rsinG)  cotGjJ  .  (12) 

Substituting  (12)  into  (9),  we  have  in  terms  of  the  following 
functions : 


y»oo 

P(t)  =  J  p(t)  p (t-x ) dt 
—  00 

(13) 

00 

P  (t  )  =-f  p (t) p (t-T  )  dt 

J  —00 

(14) 

and 

*P*( x )  =  -[  p (t)p  ( t-x )  dt 

•/ —00 

(15) 

The 

expression 

of  each  element  of  F  is  given  in  Appendix  A. 

The  lower 

bounds  on  the  variances  cf  the  estimates  for 

al' 

B,  V 

r  t  and  G  are : 

CRB  (a±)  =  (F"1)11 

(16a) 

CRB  (R)  =  (F-1)22 

(16b) 

CRB  ( x x)  =  (F“1)33 

(16c) 

CRB  (yx)  =  (f”1) 44 

(16d) 

CRB  (r)  =  (F“1)55 

( 16e) 

CRB  (G)  =  (F_1)66 

( 16f ) 
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The  results  in  this  report  are  presented  in  the  foliating  normalized  forms: 


(17a) 
(17b) 
(17c) 
( 17d) 
( 17e) 
(17f ) 


where  dx  is  the  detector  width  in  the  scan  direction.  All  results 
will  be  presented  at  SNR^IO  for  each  detector ,  where 


and 


SNR^  A  a^/oa 

4  [C  p2(t,dt] 


1/2 


=  1. 


(18a) 

(18b) 


It  is  easy  to  convert  these  results  to  any  signal-to-noise  ratio , 
SNR,  by  the  following  relationship: 


at  SNR  =  (E  at  SNR-plO)  x  (10/SNR)  (19) 

where  co  =  a,  R,  x^,  y^,  r,  and  0. 

For  the  numerical  results  presented  here,  let  us  consider  that 
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the  function  p(»)  is  generated  by  convolving  an  infinite  length 
rectangular  detector  of  width,  d  ,  along  the  scan  direction  with 
a  point  spread  function  of  the  form  [10] 

E 

sf(x,y)  =  - =■  exp 

2ttoz 

where  E  is  the  image  irradiance  at  the  center  of  the  diffraction 
pattern  and 

a  =  .3040  dx.  (21) 

The  value  of  o  is  selected  such  that  90%  of  the  total  energy  is 

collected  by  a  detector  centered  at  (x  ,y  ) . 

o  o 

Figures  6a  and  6b  present  the  Gramer-Rao  bounds  on  the 
estimation  errors  for  a  linear  array  detector  (detector  M  in  Fig. 

2)  for  the  amplitude,  a1#  and  the  in-scan  location,  x^,  respectively. 
The  upper  half  of  each  figure  is  the  Cramer-Rao  bound  computed 
for  the  1-D  problem,  in  our  previous  reports  [2]  -  [4],  where 
targets  are  located  along  the  scan  axis.  The  lower  half  of  each 
figure  shows  the  equal  performance  contours  plotted  on  the  sensor 
focal  plane  for  the  2-D  problem  with  one  target  at  the  origin  and 
the  other  at  position  (x,y).  Because  the  detector  pattern  is 
scanning  along  the  x-axis,  the  equal  performance  contours  plotted 
in  the  lower  half  of  each  figure  are  symmetric  about  the  cross¬ 
scan  axis  (y-axis).  This  syimvetric  property  is  in  common  to  all 
figures  shown  in  this  report.  Notice  that  the  equal  performance 


(x-xQ)  +  (y-yQ) 
2c^ 


(20) 
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Fig.  6a.  The  Cramer-Rao  bound  on 
amplitude  estimation  errors  for  a 
array  detector. 
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CRAMER -RAO  BOUNDS  ON  x  -FRACTION  Of  (1^ 

Fig.  6b.  The  Cramer-Rao  bound  on  the 
in-scan  location  estimation  errors  for 
a  linear  array  detector. 


contours  in  these  figures  (6a  and  6b)  are  independent  of  the  y 
component  of  the  second  target  location.  It  indicates  that  there 
is  no  way  one  can  obtain  cross' scan  position  information  from  the 
output  of  a  single  detector  alone.  In  Fig.  6b,  we  see  that  the 
CRB's  go  to  infinity  as  the  crue  separation  goes  to  zero.  This 
is  correct  for  an  unbiased  estimator;  however  (biased)  algorithms 
exist  which  produce  finite  errors  at  zero  separation.  Not 
considered  here  but  of  great  importance  for  CSO's  is  the  CSO 
detection  problem  (see  Ref.  5) .  At  small  separations,  it  is 
equally  difficult  to  recognise  that  there  are  two  unresolved 
targets  as  it  is  to  estimate  their  parameters  assuming  that  there 
are  two  targets.  The  singularity  in  estimation  error  at  zero 
separation  will  be  seen  to  influence  the  performance  bounds  for 
more  complicated  detector  patterns  later  in  this  section. 

Figures  7a  and  7b  present  the  Cramar-Rao  bounds  on  the  esti¬ 
mation  errors  for  a  chevron  array  (detectors  L  and  R  in  Fig.  3)  for 
the  amplitude,  a^  and  the  orientation  of  objects,  6,  respectively. 
The  upper  half  of  each  figure  is  a  plot  of  the  surface  formed  by 
the  Cramer-Rao  bound,  which  is  a  function  of  the  second  target 
position  (x,y) .  The  surface  is  cut-off  at  a  fixed  value.  The  edge 
of  the  plateau  forms  the  innermost  equal  value  contour  on  the  lower 
half  of  the  figure.  The  value  of  the  Cramer-Rao  bound  approaches 
infinity  as  the  (x,y)  position  approaches  (0,0).  Note  that  the 
maplitude  estimation  error  does  not  approach  infinity  except  at  the  origin 
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CRAMER -RAO  BOUNDS  ON  a^FRACTION  OF  aj) 

Fig.  7a.  The  Cramer-Rao  bound  on  the.  amplitude 
estimation  errors  for  chevron  detector  pattern. 


CRAMER-RAO  BOUNDS  ON  ©(FRACTION  OF  RADIAN) 


Pig.  7b.  The  Cramer-Rao  bound  on  the  CSO 
orientation  estimation  errors  for  a  chevron 
detector  pattern. 


(where  the  two  targets  are  colocated) ,  while  the  estimation  error 
for  the  orientation  of  the  objects  approaches  infinity  at  target 
orientations  aligned  with  either  one  of  the  detectors.  Since  the 
target  separation  is  zero  for  one  detector,  the  CRB  on  estimation 
errors  for  that  detector  is  infinite.  For  orientation  estimation, 
information  from  both  detectors  is  required  resulting  in  infinite 
error  for  an  unbiased  estimator  of  6.  As  with  a  single  detector, 
this  is  a  reflection  of  the  CSO  recognition  problem.  For  amplitude 
estimation,  it  is  sufficient  to  use  information  obtained  from  the 
detector  which  does  not  have  CSO  interference. 

Figure  8a  to  Figure  8f  present  the  Cramer-Rao  bounds  on  the 
normalized  estimation  errors  defined  in  Eqs.  (17a)  -  (17f) , 
respectively,  for  the  crow's  foot  detector  (detectors  L,  M,  and 
R  in  Fig  3) .  A  similar  format  to  the  results  shown  in  Fig.  7a 
and  7b  is  used.  One  thing  in  common  to  all  these  figures  is  that 
the  errors  increase  monotonically  toward  the  center  (0,0)  but  none 
of  these  parameters  has  errors  approaching  infinity  except  at  the 
origin.  As  seen  in  Fig.  3,  for  any  target  orientation,  there 
will  always  be  at  least  two  detectors  which  measure  a  nonzero 
separation.  Since  all  parameters  can  be  estimated  using  informa¬ 
tion  from  two  detectors,  there  are  no  singularites  as  was  the  case 
for  the  chevron  pattern. 

The  shape  of  the  contours  shown  in  Figs.  8a  -  8f  can  be  ex¬ 
plained  by  the  following  heuristic  argument:  The  estimation 
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Fig.  8a.  The  Cramer-Rao  bound  on  the 
amplitude  estimation  errors  for  a  crow's 
foot  detector  pattern. 
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ANGIE 


Fig.  8c.  The  Cramer-Rao  bound  on  the 
in-scan  location  estimation  errors  for 
a  crow's  foot  detector  pattern. 
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CRAMER -RAO  BOUNDS  ON  Y  (FRACTION  OF  DETECTOR  WIDTH) 

Fig.  8d.  The  Cramer-Rao  bound  on  the 
cross-scan  location  estimation  errors 
for  a  crow's  foot  detector  pattern. 
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ANCLE 


CRAMER -RAO  BOUNDS  ON  r  (FRACTION  OF  r) 


Fig.  8e.  The  Cramer-Rao  bound  on  the 
target  separation  estimation  errors 
for  a  crow's  foot  detector  pattern. 


C*AA«-RAO  BOUNDS  ON0IFRACTION  OF  RADIAN) 


Fiq.  8f .  The  cramer-Rao  bound  on  the 
CSO  orientation  estimation  errors  for 
a  crow's  foot  detector  pattern. 
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performance  for  individual  parameters  is  a  function  of  the  in¬ 
formation  obtained  by  each  detector  in  the  scan  direction.  For 
example,  the  effective  target  separation  measured  by  the  kth 
detector  in  the  scan  direction 


Axk  *  rfcosG  -  sinG  cotQ^] 


(22) 


changes  as  the  target  pair  orientation,  6,  changes.  Figure  3 
shows  the  detector  outputs  for  three  target  pair  orientations, 
namely,  0°,  45°  and  90°.  Table  I  shows  the  effective  target 
separation,  Axk,  measured  by  each  detector  depicted  in  Fig.  3 
for  the  three  target  pairs.  From  Fig.  6a  and  6b  we  can  approxi¬ 
mate  the  estimation  errors  as  a  function  of  Axk  as  follows: 


a2  (a)  =  C;L/(Axk)2 
a2 (Axk)  =  c2/(Axk)2. 


(23) 


The  overall  estimation  performance  for  a  parameter  w,  o(w),  can 
be  approximated  as  a  function  of  the  estimation  performance  for 
Axk  by  individual  detector  outputs,  °]c(Axk).  If  for  simplicity, 
we  ignore  the  correlation  between  r  and  0,  we  have 
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Table  I  also  shows  the  overall  estimation  performance  for  a^,  r, 
and  6  using  Eqs,  (23)  -  (25). 


TABLE  I 

APPROXIMATE  ESTIMATION  PERFORMANCE  FOR  r,  AND  0 


The  approximate  estimation  performance  for  a^,  r  and  0 
shown  in  Table  I  explains  the  trend  of  the  equal  performance 
contours  presented  in  Figs.  8a ,  8e,  and  8f.  Thus  in  Fig.  8a, 
estimation  performance  is  best  at  45°  orientation  and  worst  at 
90°  orientation.  It  is  expected  that  the  estimation  performance 
for  R  should  have  a  similar  pattern  to  that  in  Fig.  8a  for  a^. 
Combining  the  patterns  for  r  and  0  observed  in  Figs.  8e  and  8f, 
one  should  be  able  to  interpret  the  patterns  shown  in  Fig.  8c 
and  8d  on  x  and  y,  respectively. 

This  model  is  oversimplified  in  that  it  treats  the  detectors 
independently.  The  analysis  used  to  obtain  Figs.  6-8  accounts 
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for  correct  combination  of  individual  detectors  and  thus  these 
results  differ  from  those  in  Table  I. 

To  summarize  the  major  points  for  the  three  classes  of 
detector  patterns  considered  in  this  section 

(1)  The  linear  detector  can  measure  the  location  of  a  single 
target  only  in  the  in-scan  direction.  It  can  resolve  two 
CSO's  only  if  they  are  separated  in  the  in-scan  direction. 

(2)  The  chevron  detector  can  measure  the  location  of  a 
single  target  in  both  the  in-scan  and  cross-scan  directions. 

It  can  resolve  two  CSO's  if  they  are  separated  in  either 
direction  but  will  have  degraded  performance  if  their 
separation  direction  is  aligned  with  one  of  the  detectors. 

(3)  The  crow's  foot  detector  can  measure  the  location  of 

a  single  target  in  both  the  in-scan  and  cross-scan  directions. 
It  can  resolve  two  CSO's  if  they  are  separated  in  either 
direction. 

(4)  For  all  detector  types,  the  resolution  and  measurement 
performance  improves  with  increasing  signal-to-noise  ratio. 


IV .  BRICKWALL  DETECTOR  PATTERNS 


In  this  section,  the  detector  array  consists  of  M  rectangular 
detectors  with  0k  -  90°.  The  center  of  each  detector  is  not 
necessary  located  along  the  scan  axis  as  shown  in  Fig.  2,  For 
simplicity,  let  us  assume  that  the  output  of  each  detector  to  a 
unit  strength  point  source  located  at  (x,y)  has  the  same  functional 
form: 

Pk(t,  x,  y)  «  px(t-x+xk)py(yk-y)  (26) 

where  p  (•),  p  (•)  are  symmetric,  nonnegative-valued  functions 
x  y 

and  (x^,  y^)  is  the  center  of  the  kth  detector.  Let  p (t)  be  the 
autocorrelation  function  of  px(t),  i.e., 

P(t)  *  J  Px(t)px(T-t)dx.  (27) 

—  00 

Rewriting  (3)  in  terms  of  the  parameters  considered  in  Section 
II  with  the  above  assumptions,  we  have 

sk(t)  =  aiPx(t-x1+x]c)Py  (yk-yi)  +  a2Px(t-x1-rcos0+xk)py  (y^y^rsinB)  . 

(28) 

Substituting  (28)  into  (9) ,  we  have  ^  in  terms  of  p,  p,  p,  py, 
and  py.  The  expression  for  each  element  of  F  is  given  in  Appendix 
B.  The  lower  bounds  on  the  variances  of  the  estimates  for  a^,  R, 
xl'  r*  aRd  6  are  di-a9°nal  elements  of  the  inverse  of  the 
Fisher  information  matrix,  F,  as  in  Eqs„  (16a)  -  (16f)  of  Section 
III. 
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The  results  presented  In  this  section  are  values  of  the 
normalized  errors  defined  in  Eqs.  (17a)  -  ( 17 f )  in  the  last 
section.  All  results  are  presented  at  SNR^»10  of  target  Mo.  1. 
SNR-^  is  defined  in  (18a)  with 


(29) 


For  a  different  signal-to-noise  ratio,  the  same  formula  given  in 
(19)  applies.  The  same  point  spread  function  (Eqs.  (20)-(21))  is 
used  in  the  numerical  examples  presented  here. 

Figures  9a  -  9f  present  the  Cramer- Rao  bounds  on  the 
estimation  errors  defined  in  Eqs.  (17a)  -  (17f)  for  a  brickwall 
detector  pattern  shown  in  Fig.  10.  The  results  presented  here 
are  very  similar  to  those  shown  in  Figs.  8a  -  8f  for  the  crow's 
foot  pattern  in  terms  of  the  shape  of  the  equal  performance 
contours.  Notice  that  for  the  same  signal-to-noise  ratio  for 
each  detector  in  both  patterns*  the  overall  signal-to-noise  ratio 
of  the  brickwall  pattern  is  lower  than  that  of  the  crow's  foot 
pattern.  For  example,  for  a  target  located  along  the  scan  axis 
(x-axis),  the  ratio  is  about  1.5  to  1.  This  difference  should  be 
taken  into  account  when  a  comparison  of  parameter  estimation  per¬ 
formance  between  these  two  patterns  is  made.  The  amplitude  of 
each  detector  output  to  a  unit  strength  point  source  is  a  function 
of  the  cross-scan  separation  between  the  point  source  and  the 

*To  achieve  the  same  SNR  per  detector  requires  matching  the  blur 
circle  to  the  smaller  detector  size  in  the  brickwall  pattern. 
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CRAMER -RAO  BOUNDS  ON  INFRACTION  ON  lj) 

Fig.  9a.  The  Cramer-Rao  bound  on  the 
amplitude  estimation  errors  for  a 
5-detector  brickwall  pattern. 
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CRAMER -RAO  BOUNDS  ON  R  (FRACTION  OF  R) 

Fig.  9b.  The  Cramer-Rao  bound  on  the 
amplitude  ratio  estimation  errors  for 
a  5-detector  brickwall  pattern. 
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CRAMER -RAO  BOUNDS  ON  Y  (FRACTION  OF  DETECTOR  WIDTH) 


Fig.  9d.  The  Cramer-Rao  bound  on  the 
cross-scan  location  estimation  errors 
for  a  5-detector  brickwall  pattern. 
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Fig.  9f.  The  Cramer-Rao  bound  on  the 
CSO  orientation  estimation  errors  for 
a  5-detector  brickwall  pattern. 
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center  of  the  detector  as  indicated  in  Eq.  (26).  The  normalized 
peak  amplitdue  of  each  detector  output  is  plotted  as  a  function  of 
the  cross-scan  target  location  in  Pig.  10  for  the  detector  pattern 
shown  in  the  figure.  The  brickwall  pattern  differs  from  the  pre¬ 
vious  chevron  and  crow's  foot  patterns  in  using  amplitude  rather 
than  separation  estimates  from  individual  detectors  to  obtain 
information  on  target  separation  in  the  cross-scan  direction.  In 
spite  of  this  difference,  the  resulting  performance  is  quite 
similar  to  that  achieved  by  the  crow's  foot. 

The  brickwall  pattern  occupies  much  less  area  on  the  focal 
plane  than  does  the  crow's  foot  pattern.  However  to  achieve  the 
performance  presented  in  Fig.  9,  optimal  signal  processing  of  the 
outputs  of  at  least  5  detectors  was  assumed.  Simpler  suboptimal 
processing  may  be  carried  out  by  combining  adjacent  detectors  as 
shown  by  the  dotted  lines  in  Fig.  10.  In  this  way,  a  crow's  foot 
pattern  may  be  synthesized  from  the  mosaic  pattern. 
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V 


DISCUSSION  AND  SUMMARY 


In  this  report ,  we  have  derived  the  Cramer- Rao  lower  bounds 
on  the  variances  of  parameter  estimates  for  two  CSOs  separated  in 
both  the  in-scan  and  the  cross-scan  direction  as  measured  by  two 
classes  of  detector  patterns.  Numerical  results  are  presented  for 
a  particular  point  spread  function  and  a  fixed  detector  width  in  the 
scan  direction  for  both  detector  patterns.  The  trends  of  the  equi- 
performance  contours  of  each  parameter  of  interest  for  both  patterns 
are  quite  similar.  For  the  cases  considered  in  this  report,  the 
achievable  parameter  estimation  performance  of  the  5-detector 
brickwall  (mosaic)  pattern  is,  in  general,  poorer  than  the  3-de¬ 
tector  crow's  foot  patterns  for  targets  with  SNR  of  10  per  detec¬ 
tor.  This  is  because  the  overall  received  SNR  of  the  particular 
brickwall  (mosaic)  pattern  considered  is  lower  than  that  of  the 
crow's  foot  pattern.  To  achieve  the  same  performance  using  the 
mosaic  pattern  would  require  more  detector  columns  (>2)  as  well  as 
a  greater  signal  processing  load.  However,  the  brickwall  (mosaic) 
pattern  may  still  require  a  smaller  focal  plane  area  to  achieve 
the  same  performance  as  the  simple  crow's  foot  pattern. 

We  have  not  implemented  an  algorithm  to  do  2-dimensional  CSO 
detection  and  parameter  estimation.  No  simulation  performance  re¬ 
sults  have  been  obtained  to  compare  with  the  theoretical  bounds 
presented  here.  The  signal  processing  load  for  the  2-D  problem 
increases  as  a  power  of  the  number  of  detectors  in  the  pattern,  if 
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an  optimal  performance  algorithm  is  implemented.  For  a  sub- 
optimal  algorithm,  the  predicted  performance  presented  here  may 
not  be  easy  to  achieve.  This  is  an  area  which  requires  further 
investigation. 

General  trends  in  estimation  performance  were  summarized  at 
the  ends  of  Section  III  and  IV. 


APPENDIX  A 


.  THE  fisher  information  matrix  for  linear,  chevron , 

AND  CROW'S  FOOT  DETECTOR  ARRAYS 


sk(t)  =  alP(t-xk)  +  aiRp(t-Tk-Ak), 


(A.  1) 


where 


Tk  =  xrxk"yicotei 


(A. 2) 


Ak  =  rcose  -  rsin6cotek. 


(A.  3) 


Then,  the  derivatives  of  sk  are! 


3  sk(t) 


*  p(t-Tk)  +  Rp(t-Tk-Ak), 


(A. 4a) 


*  k(t) 

- - -  =  a^p (t-x  k“Ak) , 


(A. 4b) 


3  s,  (t)  ,  . 

- k -  =  -  a^lptt-Tj^)  +Rp(t-Tk~Ak)  ]  , 

3 


(A. 4c) 


*  -k(t- 


=  a^cotSj^ [p (t-Tk)  +  Rp(t-xk-Ak)  ]  , 


(A.  4d) 


3  st (t)  . 

- * -  -a,  R(cos9  -  sine  cot6k)  p  (t-xk-AK)  , 

3  r 


(A. 4e) 


•VCSi ! 


<A.4£) 


3  sk(t) 

- - -  =  a^R  r(sin0  +  cos0  cot0k)p(t-xk«Ak)  , 


where  p(t)  is  the  1st  order  derivative  of  p(t)  with  respect  to  t. 
Substituting  (A. 4a)  -  (A.4f)  into  (9)  and  rewriting  it  in  terms 
of  the  following  functions: 


CO 


p(x)  =  /  p ( t) p  (t-x ) dt , 

—  00 

(A. 5a) 

oo 

p(t)  =-  /  p(t)p(t—x) dt t 

(A.  5b) 

—  00 

and  .. 

p*(  x  )  a-  /  p(t)p(t-x)  dt/ 

—  OO 

(A. 5c) 

we  have 

M 

/  -■  (f^j)k  for  i /  j~l / 2 / .  >  / 6 

(A. 6) 

where 

(fu)k  =  ( 1+R2)  p(0)  +  2Rp(Ak) 

(A. 7a) 

(fl2)k  =  axp(Ak)  +  a^  p  ( 0 ) 

(A.  7b) 

(f15)k  =  a-^RfcosO  -  sin0cot0k)  p  (Ak) 

(A. 7c) 

(fig)k  =  -a^R  r(sin0  +  cos0cot0k) p (Ak) 

(A.7d) 

<f22>k  -  »Jo<°> 

(A. 7e) 

(f23>k  "  'al  5<V 

(A. 7f ) 
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< f 24 ^ k  *  ai2cote)c  PtAfc)  (A. 7g) 

(f33)k  »-a12(l+R2)*p*(0)-2a12Rp*(Ak)  (A.7h) 

(f34)k  «-cotek(f33)k  (A. 7i) 

(f35)k  ss-a12R(cos0  -  sin8cot0k)  [*p*(Ak)  +  R*p*(0)  ]  (A.7j) 

(f3g)k  «  a12Rr (sin6+cos6cotek)  [*p*(^k)+R*P*(0)  ]  (A. 7k) 

(f44}k  ="cot0k(f34)k  (A. 71) 

(f45)k  ="cotek(f35^k  (A*7lR) 

(f46)k  S=>'cot6k(f36)k  (A'7n) 

(f55)k  =-a12R2  (cos0-sin0cot6k)  2p*(0)  (A.7o) 


2  2  •  • 

(f5g)k  =  a1  R  r (cos0-sin0cot0k) (sin0+cos6cot0k) p (0)  (A.7p) 

(fgg^k  =~a12R2r2 (sin0+cos0cot0k) 2p (0)  (A.7q) 


and 


(f13}k 


<f14>k 


(f25}k 


(f26)k 


=  0 


where  p  and  *p*  are  the  1st  and  2nd  order  derivatives  of  p,  res¬ 
pectively.  Inverting  this  6x6  symmetric  matrix,  F,  one  should 
be  able  to  obtain  the  Cramer-Rao  lower  bound  on  the  estimates  as 
shown  in  (16a)  -  (16f)  of  Section  III. 
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APPENDIX  B 


THE  FISHER  INFORMATION  MATRIX  FOR  BRICKWALL  DETECTOR 
PATTERNS 


Let 


sk(t)  =  aiPx(t-Tk)py(yk“yl)  +  alRpx(t“Vrcos*)py(yk“yl“rsin<,,) 


(B.l) 


where 


Tk  -  w 


(B.  2) 


Then,  the  derivatives  of  sk  ares 


3  sk(t) 
3  a. 


-  Px(t“Tk)py<yk“yl)  +  RPX(t"Tk-roos0)py(yk-yl“rsine)  ' 


(B. 3a) 


3  sk(t) 


3  R 


=  a1Px(t-  yrcos0)  p  (y^-y^-sine)  , 


(B. 3b) 


3  sk(t) 
3  x. 


■  (yk-y1)-»1*Px(t-Vra»e)py(y|t-y1-r.ine) 


(B.3C) 


3  sk(t) 
3  y1 


=  -aipx(t-xk)p  (yk-y1)-a1Rpx(t-Tk-roos0)py(yk-y1-rsine)  , 


(B. 3d) 


3  sk(t) 


=  -a^Rcos0px (t-Tk-rcos0) p^ (yk-y^-rsin6) 

-  a1Rrsin0px (t-Tk-rcos0) p^ (yk-yj-rsin0)  , 


3  r 


3  s^(t) 


3  0 


“a^Rrsin©^  (t-ik-rcos0)py  (y^-S^-rainO) 

-a1Rrcos0px(t-Tk-rcos6^py (y^-y^-rsinG) $ 

(B.3f) 


where  p  (•)  and  p  (•)  are  the  1st  order  derivatives  of  p  and  p  , 
x  y  x  y 

respectively.  For  notational  simplicity,  we  will  use  p^  and  P2 

in  the  rest  of  this  Appendix  to  denote  p  (y,-y, )  and  p  (y  -y  -rsin9) , 

y  k  x  y  k  1 

respectively. 

Substituting  (B.3a)  -  (B.3f)  into  (9)  and  rewriting  it  in 
terms  of  p,  the  autocorrelation  function  of  px  defined  by  (11), 
and  its  derivatives  p,  and  *p\  we  have 


M 

if  Z  <*„>„. 

No  k=l  13  K 


for  1, j  *  1,2,.. .,6 


(B.  4) 


where 

(fi;i)k  =  (p12+p22R2)  p  (0)  +  2p1P2Rp  (,rcos6) 

2 

*f12*k  =  alplp2p (rcos0)  +  a1RP2  P<°) 

(fl4Jk  =-al(plpl+R2p2p2) p(0)  ”  a1R(p1P2+p1p2) p (rcosO) 

(fl5) k  =  aiRpip2  cos9P(rcos0)  “  a^Rp^PjSinOp (rcos6) 

-a1R2p2p2sin0p(O) 

(f 16> k=”alRrplP2  sin0P(rcose)  ”  a1Rrp1p2cos0p (rcos0) 

-a1R2rp2p2cos0p (0) 


(B.5a) 
(B. 5b) 
(B. 5c) 

(B.5d) 

(B.5e) 
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(B.5f) 


<Vk  ‘  ,llS2p(0) 

(f23}k  ““ a12plp2^(rCOB°) 

(f24)k  »-a12p1P2P(rcose)  -  a12Rp2P2P(0) 

(f25) k  ““ ■a12*PABiBep(0) 

(f26)k  =-a12Rrp2P2co80p(O) 

(f33)k  =-a12(p12+P22R2)P(0)  *  2a12Rp1P2V(rcoa8) 

(f34)k  -  «i2r(PiP2  ’  PiP2)P(rcose) 

(f35)k  »-a12Rp1P2cose*p’(rcos9)  -  a1‘R2p22cose*p  (0) 

+  a12Rp1P2sin0p(rcos6) 

(t36)k  =  a12Rrp^P2sin6p*(rcos6)  +  a^Rrp^ooseS  (rcos0) 

+  a12R2rp22sin6p*(0) 

(*44)k-  (al2p12+a12R2p22)P(0)  +  2a22Rp^P2P (rcos0) 

(tt5)k  =.-a12RplP2oos6p(rcose)  +  a12Rp1P2sin6p  (rcos0) 

+  a12R2p22sinBp (0) 

(£46)k  =  a12Rri1P2aineJ(rcos0)  +  a^Rrp^eosep  (rcose) 

+  a12R2rp22cos8p  (0) 

(f  55)  k  -“ai2R2  [cos28p22*p(0)  "  sin2eP2  p(0)] 


(B.5g) 
(B.  5h) 
(B. 5i) 
(B. 5 j ) 
(B.  5k) 
(B. 51) 

(B.  5m) 

(B.5n) 

(B.5o) 

(B.  5p) 

(B.5q) 
(B. 5r) 


48 


t  »  ; 


(f56}k 


a12R2r8in8cos^P22*p^0^  + 


(f66)k  — a12R2r2^p22sin26*p*(0)  -  p22cos26p (0)J 


(B.5s) 


(B. 5t) 


.1 
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